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I . INTRODUCTION 


The nature of the local atmosphere of lo and the manner in which 
this atmosphere escapes to produce extended gas clouds and to supply 
the heavy ion plasma torus threading the satellite orbit are the central 
topics of this research, lo is known to have extended gas clouds of 
neutral sodium, potassium and oxygen detected by ground-based observa- 
tions. In addition, detection of a neutral sulfur cloud has very 
recently been reported based upon UV emission data obtained from a 
rocket flight (Durrance, Feldman and Weaver, 1983). lo is also known 
to have a dense hot plasma torus composed primarily of oxygen and sul- 
fur ions discovered by ground-based observers and more recently measured 
by Instruments of the Voyager spacecrafts. The two basic goals of this 
research are to characterize the satellite emission conditions operative 
at lo for the sodium, oxygen and sulfur neutral gas clouds and to help 
characterize the satellite-ion source and magnetic diffusion of ions in 
the near lo environment. 

To achieve the two research objectives, two different approaches 
are being followed. The first is to identify the atom escape character- 
istics of sodium from the satellite by analysis of the substantial neu- 
tral sodium cloud data base that has been obtained from Earth-telescope 
observations over the past ten years. The second is to explore by model- 
ing analysis the implications of the rather recently discovered lo oxygen 
cloud (Brown, 1981) and the very recently reported detection of an lo 
sulfur cloud (iDurrance, Feldman and Weaver, 1983). The discovery of the 
oxygen cloud in 1981 provided the first direct link between the loss of 
local atmospheric gases form the satellite and the supply of 0^ ions to 
the plasma torus. The relative importance of direct ion escape (i.e., 

S , 0 , SO , SO 2 , etc.) and ionization of neutral gas clouds (i.e., S, 0, 

SO,, etc.) as ion supply channels for the plasma torus is currently not 

^ + + ' 
known. Better definition of the net rate at which 0 ions and S ions 

are supplied by ionization of the neutral gas clouds has been one of the 

major thrusts of our research effort this .year. In addition to this net 

plasma source, the neutral gas clouds have also recently been realized 

(Johnson and Strobel, 1982) to directly alter to a significant extent the 

relative abundance of the plasma torus ions through neutral-ion exchange 
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reactions. The spatial distributions of the neutral clouds are in turn 
controlled by loss processes involving collisional ionization and charge 
exchange reactions. This complex coupling between the lo plasma torus 
and the neutral clouds has been a subject of Intense Interest in our 
research effort this year. 

Progress made toward our two research objectives in the second year 
of this endeavor is reported here. The analysis of the sodium data has 
received less emphasis than the exploratory modeling of the oxygen cloud. 
Significant advancement has however been made toward each objective and 
is discussed in the following section. 



II. PROGRESS DURING THE SECOND YEAR 


Introduction 

The strategy adopted during the second year has been to focus 
more effort upon the exploratory modeling of the lo atomic oxygen and 
atomic sulfur clouds and less effort upon the analysis of the lo so- 
dium cloud data. This strategy was adopted to optimize our scientific 
program in response to reduced budgetary support. The sodium data 
analysis effort has thus been restricted to acquisition of new data 
and to determination of suitable modeling basis functions. Quantita- 
tive analysis of the sodium data will be initiated in the third year . 

The three-year plan for our data analysis modeling is summarized in 
Table 1. 

Modeling of the lo Oxygen Cloud 
Model Calculations 

Further exploratory model calculations for the oxygen cloud, in- 
corporating Improvements in the lo plasma torus electron data adopted 
during the first year, have been performed in the second year. The 
distribution of 01 in the extended cloud was placed on an absolute 

O 

scale by fixing the modeled 6300A emission intensity to the observed 

O 

value (Brown, 1981) , and model predicted intensities for the 01 1304A 

O 

and 880A emissions were below observed upper limit values (Moos and 
Clarke, 1981; Shemansky, 1981). The loss process controlling the 
neutral population in these calculations was ionization by lo plasma 
torus electrons, and model results appropriate to the plasma conditions 
at the time of the Voyager 2 encounter with Jupiter are summarized in 
the second column of Table 2 (for more details see Smyth and Shemansky, 
1983 Included in the Appendix) . 

The effects of Including a neutral sulfur cloud and also of includ- 
ing charge exchange reactions between plasma torus ions and the neutral 
01 and SI cloud atoms were evaluated and found to be significant. A 
rough estimate to Include these effects was made by scaling (using as 
a guide the plasma torus calculations of Brown, Shemansky and Johnson, 
1983) our lo oxygen cloud results (listed in the second column of Table 
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Table 1 

Three-Year Plan for Data Analysis Modeling 
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SUMMARY OF MODEL RESULTS 
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2) and assuming an lo sulfur flux of one-half the scaled oxygen flux. 
These scaled model results are summarized In the third column of Table 
2 and are discussed in more detail by Smyth and Shemansky (1983) . The 
Ion loading, plasma mass loading and Ion energy Input rates, estimated 
In Table 2, strongly Influence the composition, diffusion, field-aligned 
current systems and energy budget of the plasma torus and are of funda- 
mental Importance to models of the Jovian magnetosphere. The cloud mass 
loss rate, estimated in Table 2, is of fundamental importance in under- 
standing the local atmosphere of lo and its interactions with the plane- 
tary magnetosphere. Obtaining a better estimate of these rates and 
their spatial character is therefore very important. 

Model Improvements 

Two major model Improvements not Included In the model results 
discussed above, were Initiated in the latter part of the second year: 

1. Introduction of the oscillating motion of the plasma 
torus about the satellite plane, and 

2. explicit inclusion of charge exchange reactions in 
the model, which affect both the cloud atom life- 
time and the model predicted values of the physical , 
rates listed in Table 2. 

The first improvement has been implemented and most of the work required 
to implement the second improvement has been undertaken. To complete 
the second improvement, the spatial density of ions, in the torus must 
be specified in two-dimensions. This information is at present only 
partially available. Parallel efforts have however been initiated to 
develop the necessary description of the ion densities in the plasma 
torus from the diverse and sometimes apparently contradictory sets of 
Earth-based, rocket, lUE, and Voyager measurements. Recent analysis 
of Voyager ion data by Bagenal (1983) has provided the required ion 
density information for the radial interval from 4.9 to 5.4 Jupiter , 
radii. Additional information for the ion densities for larger radial 
displacements is also being sought from the analysis of Voyager UVS 
' data by Shemansky (1983) . The UVS inferred ion densities will be 
used to complement the Voyager ion data which does not uniquely define 
the ion densities in the hot torus (>5.7 Jupiter radii). 
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Charge exchange reaction rates for the model have been obtained 
from Johnson and Strobel (1982) and Johnson (1983) . The set of rele- 
vant reactions adopted In the model are listed in Table 3 and the 
velocity dependence of the charge-exchange reaction rates have been 
Included in the model as prescribed by Johnson (1983) . Relevant 
electron Impact ionization and recombination reactions are listed in 
Table 3 for the neutral species. Ion-ion and electron-ion reactions 
not producing neutral species as products are excluded from Table 3 
since they are not required to calculate the model results in Table 2. 
Additional model results, such as the production or loss rate of a 
particular ion species, would however require these additional reac- 
tions to be incorporated. New model results will, in the near future, > 
be limited to calculating improved values for only the rates listed 
in Table 2 (see the Appendix for further discussion) . 

Modeling of the lo Sulfur Cloud 
Model Calculations 

The lo sulfur cloud model, developed during the first year of the 
data analysis was used to calculate the intensity on the sky plane of 

O 

electron impact excited emission lines at wavelengths of 11,306A, 10, 
82 oA, 1115 %. and 4589A. Model results for the stronger optical emission 

O 

line at 7725A are shown in Figure 1 for the same observational geometry 

O 

used in the 6300A oxygen model calculation (simulating the detection 
measurement of Brown in Figure 4 of the Appendix). The 1125 % intensity 
contours, in Rayleighs, are about an order of magnitude dimmer near lo 

O 

than the 6300A intensity contours for the oxygen cloud and are fixed 
by assuming a sulfur escape flux of one-half the value determined for 
the oxygen cloud. The sulfur cloud is more tightly confined near lo 
because its electron impact ionization lifetime shown in Figure 2 is 
about 4 to 5 times smaller than the electron impact ionization life- 
time of oxygen (see Figure 3 in the Appendix) . This is illustrated 
more clearly by comparison of the column density of both clouds as 
viewed from above the satellite plane in Figure 3. 

O 

Near lo, the intensity of the other optical line at 4589A is 
about five times dimmer than the 1125 % line and therefore will not be 
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Table 3. Charge-Exchange and Electron Impact Reactions for 
Neutral Oxygen and Sulfur 


Reaction 


1. o'*' + 0 0 + o'*' 

2. o'*' + S -V 0(^1>) + s'*'(^P) 

-V 0(^D) + S'*'(^D) 

3. s'*" + S -> S + s'*" 

4. s'*' + 0 S + o'*" 

5. s'*^ + S S + S'*^ 

+ 4 + 4 4 

^ S ( S) + S (3s, 3p ; P) 

6. S"*^ + 0 -»■ S‘*'(^P) + Q'*'(^D) 

+ 2 +2 

S ( D) + O^(^D) 

+ 2 +2 

S ( D) + O^(^P) 

+ 4 +2 

S ( S) + 0 ( P) 

7 . 0“^ + 0 -»■ 0 + O"*^ 

0'*'(^S) + 0'*'(2s, 2p^; S) 

0 (‘‘d) + 0^(2s, 2p^; . P) 

8. O'*^ + S 0'*'(^S) + (S'*')* 

^ O-^(^D) + (S-^)* 

^ O^(^P) + (S^) 

^ 0(^P° 3d, ^P®) + s"*^ o'*' + s'*^ + e 

9. S'*'*^ + 0 ^ s'*'*'(3d, ^P) + o'*'(^S) 

-► S’*^(4s, ^P) + 0'*'(^S) 

S'*'(^P) + 0'*^(^P) 

S'*’(^D) + 0'*'*'(^D) 

-> S‘*'(''S) + 0-*^(^D) 

10. s'*'^ + S -)- s'*'*'(4p, ^P) + s'*'(^S) 

-)■ S'’^(3d, ^P) + S'*'(^P) 

(S"*^)* + S'*'(^D) 

(s-^)* + S'*^ 

11. e + 0 -i- o'*” + 2e 

12. e + S s'*" + 2e 

13 . o'*" + e 0 + hv 

14. s"*^ + e -»■ S + hv 
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Vertical Displacement from the Centrifugal Equator (Rj) 


SI Electron Impact Ionization Lifetime in thelo Plasma Torus 

(Voyager 2 Plasma Conditions) 



Radial Displacement from Jupiter (Rj) , 

Figure 2 . The lifetime is calculated using the 
plasma conditions adopted for the 
Voyager 2 encounter with Jupiter. 
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lo; Atomic Oxygen Torus lo: Atomic Sulfur Torus 
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considered further. The IR line at 10,820A is, however, comparable to 
the 6300A oxygen line emission intensity near lo while the IR line at 

O 

11,306A is about a factor of 3.5 dimmer. Comparable or brighter emis- 
sion lines at UV wavelengths of 1304A, 1425A, 1485A and 1820A, also 
exist for neutral sulfur, as Illustrated in Figure 4 for only one 

Q 

temperature of 80,000“ K. The IR emission at 10,820A provides the best 
line for ground-based detection of neutral sulfur while the UV emission 
at 182 oA provides the best line for rocket or Earth-orbiting satellite 
detection. Appropriate UV and IR observers have been notified of these 
emission lines and interest in observing the neutral sulfur cloud is 
growing . 

Model Improvements 

The model Improvements discussed earlier for the oxygen cloud model 
have also been Implemented for the sulfur cloud model. Progress paral- 
lels that achieved for the oxygen model. In addition, arrangements 
have been made with Shemansky (1983) to obtain approximate expressions 
for the electron emission rates for the UV lines of atomic sulfur listed 
above. Computations of the intensity of these UV emissions on the 
sky plane will be performed during the third year. 

Analysis of the lo Sodium Cloud Data 

The quantitative analysis of the lo sodium cloud data has been 
divided into five stages of activities which are summarized in Table 
4. Due to the reduced budgetary support available, the original pro- 
posed scope of this data analysis modeling effort cannot be performed. 

A data set of more limited scope has therefore been selected with the 
overall goal of our modeling effort remaining unchanged. In order to 
retain the virtue of the analysis in step 4 and step 5, the limited 
data set will consist of both spectral data (the sodium line profile 
data) and some spatial intensity data (sodium D-llne intensity varia-' 
tions on the sky plane as Imaged through an observing slit). In the 
first year, efforts were directed to stage 1 activities of Table 1. 
Available line profile data [Trafton, 1975; Trauger, Roesler and MUnch, 
1976; Trauger, 1977; Trafton and Macy, 1977; Trafton, 1981] and spa- 
tial intensity data [Brown et al., 1975; Trafton and Macy, 1975; 
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ELECTRON IMPACT EXCITATION RATE (cm'* sec') 


ELECTRON TEMPERATURE (ev) 

0.1 1 10 100 



ELECTRON TEMPERATURE (»K) 


Figure 4. These electron Impact excitation rate coefficients for 
neutral oxygen and sulfur were computed by Shemansky 
(1982) . The visible and IR emission rates were based 
upon close coupling calculation taking Into account 
configuration Interactions that were performed by 
Henry (1981). 
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Trafton, 1977; Brown, 1977; Trafton and Macy, 1978; Mekler, 1982] were 
acquired and preliminarily evaluated. In the second year, efforts were 
directed to stage 2, with priority given to the line profile data. 

The lo sodium cloud model, improved to include the space and time 
varying electron impact ionization sink of the oscillating plasma 
torus, will be used to generate the basis functions for analysis of 
the data set in the third year (steps 3-5) . The purpose of the second- 
year effort (stage 2) has been to specify the physical model parameters 
(i.e., the exospheric escape conditions of sodium and the plasma torus 
sink conditions) sufficiently accurately so that the basis functions 
computed by the model are physically appropriate for data inversion. 
Improvements in the description of the electron density and temperature 
have been Implemented in the model this year. Progress in selecting 
the exospheric escape conditions has resulted from (1) studying the 
character of the line profile data independent of the model and (2) 
applying new model results that have very recently become available 
from research performed in one of our other NASA supported projects 
(NASW-3387). The new model results and their implications for this 
project are briefly discussed below. 

In our NASA supported research project (NASW-3387), modeling of 
the sodium directional features discovered by Pilcher (Hartline, 1980) 
has been in progress for the last three years. This modeling has 
recently resulted in the discovery of new and direct evidence for a 
magnetospherlc-drlven gas escape mechanism for lo. , This modeling 
analysis has used the lo sodium cloud model with the oscillating 
plasma torus. It has shown that the directional features outside of 
lo's orbit can be understood in terms of (1) an oscillating plasma 
torus sink, (2) a high-velocity (•'20 km sec ^) sodium source whose 
direction must be nearly perpendicular to lo's orbital motion, and 
(3) a non-uniform flux of sodium from lo that is enhanced near the 
equator of the satellite. Adopting these exospheric escape conditions, 
preliminary efforts indicate that such high-velocity sodium atoms may 

t 

also have the correct characteristics to explain the asymmetric skirt 
exhibited in the sodium line profile data. Basis functions consistent 
with these escape conditions will be evaluated more fully in the third 
year and may prove to eliminate past ambiguities that have existed in 
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choosing unique basis functions for data inversion of the line profile 
data. These escape conditions are also Important in properly under- 
standing the north-south alternating spatial intensity data [Trafton 
and Macy, 1975; Trafton, 1977] which comprises part of the complemen- 
tary (i.e., non-spectral) information in our limited data set. 
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Abstract 

The escape of atomic oxygen from lo and its distribution in the 

magnetosphere of Jupiter is described in a model calculation. The loss 

process controlling the neutral population in this calculation is ionization 

by plasma electrons. The 01 distribution is placed on an absolute scale by 

fixing the modeled 6300 % emission intensity to the observed value. Predicted 

0 0 

intensities for the 01 1304 A, and 880 A emissions are included in the 

calculation. The model indicates a satellite emission flux of 

1.5 X 10* cm-* s-*t an ion loading rate of 6 .2 x 10** ions s-*, an oxygen ion 

mass loading rate of 16.6 kg s-*» and an oxygen ion energy input rate of 

2.7 X 10** watts. The effect of including a neutral sulfur cloud and charge 

exchange reactions in the plasma torus is significant. A rough estimate 

including these effects raises the source flux for oxygen atoms to 

1 .2 X 10** cm-* S-*. The inclusion of an assumed satellite flux for sulfur 

atoms of 6 .0 X 10* cm-* s-* then produces an ion loading rate of 

~4.0 X 10*^ ions s-*> an ion diffusive loss time of ~200 days, a plasma mass 

loading rate of ~150 kg s-*, a satellite mass loss rate of ~270 kg s-*, and a 

maximun ion energy input of ~4 x 10** watts. The energy lost by radiation 

(3 X 10** yn) from the Voyager 2 epoch plasma torus exceeds the estimated 

energy input through the production of now ions by an order of magnitude. 



1 . Introduction 


lo, the innermost Galilean satellite of Jupiter, is known to be the 
source of a substantial plasma torus in the Jovian magnetosphere composed 
primarily of oxygen and sulfur ions. The satellite is also known to have 
neutral gas clouds extending spatially beyond its local gravitational control. 
To date, neutral gas clouds of sodium (Brown 1974), potassium (Trafton 1975; 
Munch, Trauger, and Roesler 1976; Trauger, Rocsler, and Munch 1976) and atomic 
oxygen (Brown 1981) have been detected from ground-based facilities, and at 
least one other species, neutral sulfur, has now been reported as an observed 
emission by Durrance, Feldman, and Weaver (1982) . 

The role of the neutral gas clouds of lo as controlling factors for the 
maintenance and morphology of magnetospheric plasma has not been generally 
recognized. However, if mass is continuously lost from the plasma as it must 
bo through various processes, neutral particles must ultimately serve as 
fodder for the maintenance of the system. If neutral clouds exist in the 
space of the plasma torus they must certainly be a source of .new ions, The 

t 

recent discovery of the atomic oxygen cloud in the plasma torus remote from lo 
(Brown 1981) demonstrates the presence of a . significant source of both plasma 
and energetic neutral atoms for the magnetosphere because of the recent 
realization that neutral-ion charge exchange reactions in addition to electron 
collisions^are controlling factors in the lo plasma torus (Brown, Pilcher and 
Strobcl 1983; Johnson and Strobel 1982; Drown, Shemausky, and Johnson 1983; 
Brown and Shemansky 1982) . The action of these processes is such that the 
characteristics of the lo atomic oxygen cloud and the properties of tl^e 
magnetosphere are strongly coupled. 

The significance of the recent detection by Brown (1981) of an atomic 
cloud near the orbit of lo, achieved by observing its electron-excited 
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emission in the 6300 A line, is that it provides ns with the first direct 
information to link the escape of neutral gas from the satellite and the 
observed oxygen and sulfur ions in the plasma torus. The population of such 
neutral gas clouds in the Jovian circumplanetary space is governed by the 
source rate at lo and by the subsequent loss processes which occur for the 
escaping atoms along their orbital paths in the magnetosphere. The 
interaction of these clouds with the plasma torus forms the dominant sink, 
converting the atoms to ions in collision with electrons and existing ions. 

The observation of excited emission from the 01 atoms through their co^ision 
with the torus plasma electrons can then provide a means of estimating both 
the neutral cloud density and the resulting mass, ion and energy input 
supplied to the magnetosphere by the cloud. Those quantities, which are 
fundamental to understanding the connection of the plasma torus to the 
magnetosphere, haye been the subject of heated discussion in the literature. 
The torus mass loading and diffusive loss rates have been estimated by a 
variety of methods (Sullivan and Siscoe, 1981; Eviatar and Siscoe, 1980; 
Richardson £t £1. , 1980; Richardson and Siscoe, 1981; Dossier, 1980; Brown, 
1981; Shemansky, 1980; Thorne, 1981; Hill, 1980) providing values separated by 
factors of 100 or more. An uncritical reader of the various arguments 
therefore cannot arrive at a sensible conclusion in regard to questions 
relating to how energy is delivered to the torus for its maintenance, and the 
nature of the interaction of the torus with Jupiter's atmosphere and 
magnetosphere. , 

The emphasis of this article will be to explore what can be learned about 
the escape processes of atomic oxygen from lo and the plasma source that the 
neutral gas forms in the planetary magnetosphere. Other plasma sources for 
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the nagnetosphere such as direct ion escape from Io« escape of molecules such 
as SO, followed by dissociation and ionization, or escape of molecular ions 
such as S 02 ^ followed by dissociative processes are also possible, but none of 
these processes have yet been positively identified. To explore the 
characteristics of the atomic oxygen cloud and its impact on the 
magnetosphere, a model for the gas cloud has been developed and suitable 
calculations performed. 

The model calculation requires establishment of a reasonably accurate 
temperature and density structure for the plasma, which we obtain from a 

I 

combination of Voyager i£ situ and EUV data. The loss process providing the 
sink for neutrals in this calculation is electron impact ionization only, 
although we realize that charge exchange reactions are an important additional 
loss process in some regions of the torus, however, pending further detailed 
calculations, we make rough estimates below of the effect of ion-atom 
reactions on the atomic oxygen distribution and source rates., based on recent 
work by Brown, Shemansky, and Johnson (1983) . 

Qualitatively the Brown (1981) observation of 01 6300 X emission 
indicates a rather low abundance implying low mass loading and source rates 
compared to most of the earlier published results. The energy injection rate 
based on the calculated production rate of ions is too low to compensate for 
the known radiative cooling rate. This result is uot new in the published 
literature since Shemansky (1980) and Shemansky and Sandel (1982) have 
obtained similar results based on other considerations. Bowever therje are 
conflicts with the other work indicated above, and the present results are by 
no means generally accepted. We present the model calculation below and cite 
the arguments and uncertainties in the issue. 
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2 . Atomic Oxygen Cloud Model 

Brovirn (1981) by adding six individual observed spectra obtained an 
intensity of 8+4 Rayleighs in the 6300 A emission lines. This difficult 
detection, which required three hours of integration time, was made possible 
by the high spectral resolution of his instrument which separated the lo 
emission line from the bright telluric oxygen airglow line of approximately 80 
Rayleighs. The lo related emission in the 6300 X lines is excited by 
inelastic impact of oxygen atoms and lo plasma torus electrons (Brown. 1981) . 
Other attempts to observe electron-impact excited emissions of the oxygen 
cloud at wavelengths of 1304 X and 880 A. summarized in Table 1 have provided 
only upper limits. The plasma torus in addition provides a sink for the lo 
atomic oxygen cloud through electron impact ionization and ion charge exchange 
reactions. The effect of charge exchange processes, although important in the 
dense regions of the hot torus near Io*s orbit (Brown, Shemansky and Johnson. 
1983) and in the inner cool region where the torus electrons become too cold 
for impact ionization (Johnson and Strobel. 1982) has not been included here 
as mentioned earlier. The brightness of the oxygen atom emission in a given 
volume element of the cloud is thus determined by the local electron number 
density and temperature, while the number of oxygen atoms in that volume 
element is determined by the previous ionization time-history of these atoms 
in the planetary environment and the source injection rate. 

and 880 X emission intensities, and 
instantaneous oxygen ion creation rate of the oxygen cloud is achieved by 
following in circumplanetary space the trajectories of many oxygen atoms 
ejected from lo and by determining proper weights along these trajectories to 
reflect their real time volume excitation and ionizationrates. The 


Modeling the density. 6300 X. 1304 X. 
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trajectories i& circiimplanetary space from the assumed isotropic ejection 

I 

process at lo are determined by solving the circular restricted three-body 

I 

equations of motion for each atom with suitably specified initial conditions. 
This procedure has been documented in detail by Smyth and McElroy (1977 and 
1978) who applied this method to modeling the lo sodium cloud. To determine 
the proper excitation and ionization weights along these trajectories, the 
volume excitation and ionization rates of atomic oxygen must be expressed as a 
function of the plasma parameters, and the plasma parameters must also be 

specified as a function of spatial coordinates. 

) 

The description of the plasma torus that is adopted for modeling purposes 
in this paper will be limited to the electrons. Ion-ion and ipn-neutral 
collision processes will be included in future modeling when improved spatial 
density information for the ion species becomes available. We adopt an 
electron density-temperature spatial distribution based on Voyager 
experimental data. A single Blaxwollian electron energy distribution is 

I 

applied although we realize such a description is not strictly correct since 
it is known that the' electrons arc uon-Maxwellian and that the microscopic 
state varies with radius and magnetic latitude within the torus (Scudder. 
Sittler and Bridge. 1981) . This non-Maxwell ian distribution is actually 
characterized by a dominant cool component (S-26ev) and a small hot component 
(626-1200 ev) . The relative abundance of the hot component measured at 
Voyager 1 encounter increases with radius, having values of 0.02%. 1.38% and 
7.93% for radial values of 5 , 5 , 7.8 and 8.9 Jupiter radii (R.) respectively 
(Scudder, Sittler, and Bridge. 1981) . Further, observations of an upper limit 
of OIII abundance in the hot torus (Brown. Shemansky and Johnson, 1983) have 
placed a stringent upper limit on the high energy electron component in the 
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post Voyager 2 encounter epoch. According to the Brown. Shemansky and Johnson 
calculations, which include charge exchange reactions in the central region of 
the torus, the limitation on the population of hot electrons is relatively 
independent of ion diffusive loss time. Quantitatively, a relative hot 
electron abundance of only 0.02% is above the limit indicated in the latter 

work. On this basis, the contribution of the hot component to the electron 

' 0 0 0 
Impact ionization and to the 6300 A. 1304 A and 880 A excitations of atomic 

oxygen is not significant compared to the cool component contribution. 

Although ionization rate coefficients are much higher for the hot electron 

component, the density relative to the cool component is too low to contribute 

0 

significantly to the lifetime (see Table 1) . Excitation of 6300 A emission is 

0 

inefficient at high temperatures (Table 1) whereas excitation of 1304 A and 
880 X is higher (Table 1) but is offset by the low density of the hot 
component. Because of this, specification of only the low temperature 
electron component of the plasma is required in the modeling' work to follow. 

Figure 1 shows the two radial temperature profiles that have been adopted 
in this paper for the cool electrons in the torus appropriate to the two 
encounter conditions of Voyager 1 and Voyager 2 with Jupiter. The ion 
temperature profile for the Voyager 1 encounter condition is also indicated 
for comparison, as deduced from the in situ direct plasma measurements by 
Bagenal and Sullivan (1981) using their common temperature model. For radial 


displacements from Jupiter less than 4.9 R.. the adopted electron temperatures 

J 

in Figure 1 are assumed to have a constant value of 5.8 z 10’ K (O.S ev) . 


whereas for radial displacements from 4.9 R. to 5.3 R. values equal to the 

J J 

measured ion temperature are assumed. For radial distances from 5.3 R^ to 

7.0 R.. the electron temperatures are equal to average electron temperatures 
1 
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deduced from analysis of the extreme ultraviolet emission from the torus ions 

obtained by the UVS instrument of tlie Voyager 1 and Voyager 2 spacecraft. In 

the 7.0 to 8.5 Rj region, temperatures are interpolated using the in situ 

measurements of Scudder, Sittler and Bridge (1981) for the Voyager 1 epoch, at 

7.8 and 8.9 Rj. Scudder, Sittler and Bridge (1981) also provide a measured 

temperature at S.S Rj which lies approximately 20 percent below the 

interpolated curve in Figure 1. Beyond 8.5 R. the electron temperatures are 

J 

assumed to have a constant value of 3 x 10* K (25.8 ev) . The exact value of 
the temperature beyond 8.5 R. is not important since ionization rates show 
little dependence in this temperature region. The Voyager 2 epoch curve 
beyond 7.0 Rj is established by assuming the terminal temperature at 8.5 Rj 
was the same as the measured value at the Voyager 1 encounter. 

The radial profile of the electron number densities in the centrifugal 
equator of the lo plasma torus, appropriate to the encounter conditions of 
Voyager 1 and Voyager 2 with Jupiter, is shown in Figure 2. The two electron 
density profiles from the Voyager 1 ^ situ plasma data (Bridge, Sullivan and 
Bagenal, 1980; Bagenal and Sullivan, 1981; Bagenal, 1981), deduced by equating 
them to the magnitude of the total ion charge density for the common 
temperature and thermal speed models, are shown between the radial distances 
of 4.0 R. and 8.5 R. . The two profiles are identical for radial displacements 

J J 

less than 5.8 R., and an extrapolated value has been assumed between 4.0 R. 

J J 

and 4.9 R. based upon refined data at 4.9 R. and 5.0 R. provided by Bagenal 
1 ' 11 


(1981). The average electron number density profiles obtained from analysis 
of the ultraviolet ion emission data from the plasma torus measured by the UVS 
instruments of Voyager 1 and Voyager 2 are also indicated and have been 
enhanced by a multiplicative factor of 1.25. This enhanpement factor has been 
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chosen to adjust the spatially averaged values of the Voyager 1 electron 
density determined from the UVS data to their maximum values in the 
centrifugal plane as defined by the Voyager 1 in situ plasma data. For model 
calculations appropriate to the Voyager 1 encounter conditions, the adopted 
number densities in the centrifugal plane are defined by the ^ situ results 
from 4.0 Rj to S .8 Rj and by the enhanced UVS results for radial displacements 
larger than S.8 Rj. For model calculations appropriate to the Voyager 2 
encounter conditions, the adopted number density in the centrifugal plane will 
be defined by the in situ Voyager 1 results increased by a multiplicative 
factor of 1.4 in the radial interval from from 4.0 R^ to 5.8 Rj and by the 

modified UVS results of Voyager 2 for radial displacements larger than 5.8 R.. 

J 

In the model calculations described below, the electron number density in 
the lo plasma torus is specified in two-dimensions, n(r.z). having radial (r) 
and vertical (z) dependencies above and below the centrifugal plane. The 
number density is assumed to be independent of oiagnetic longitude. The • 
two-dimensional electron number density profiles in the centrifugal plane are 
determined from the radial number density profiles in the centrifugal equator. 
n(r) adopted in Figure 2, by an approximate form of a scaling law utilized 
earlier by Bagenal and Sullivan (1980). 

n(r.z) = u(r) , (1) 

In this expression the distance along the magnetic field line from the 
centrifugal plane has been approximated by the vertical distance from the 
centrifugal plane. This approximation is acceptable for the small vertical 
departures of interest in the lo plasma torus. The scale height H. in units 

of km. is defined by 

z T 

H «= (1.7924 X 10») 

M 


( 2 ) 
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where the tcriperature T, in units of has been chosen as the cool component 

electron temperature profile adopted in Figure 1. The value of the effective 

Mass M, in units of A!dU, has been selected to bo 8.2 by requiring the 

calculated number density n(r,z) for the Voyager 1 encounter conditions to be 

ip reasonably good agreement with the more carefully described twp-dimensional 

electron number density distribution reported by Bagcnal and Sullivan (1981) . 

having specified the spatial dependence of the electron temperature and 

number density in the torus, the electron impact ionization lifetime add the 
0 0 0 

6300 A, 1304 A and 880 A volume excitation rates for atomic oxygen can be 
determined using the data in Table 2. The excitation rates given in Table 2 
were calculated using the Brook. Harrison and Smith (1978) measurements of the 
01 ionization cross-section and an 01 electron excitation model calculated by 

Shemansky (unpublished) based on available measured and theoretical collision 

0 

strengths (see Table 2) . The 6300 A collision strengths, in particular, were 
based upon the Henry. Burke and Sinfailom (1969) calculations. Equivalent 
analytic expressions were used in the model calculations. The two-dimensional 
electron-impact ionization lifetime for oxygen in the torus is illustrated in 
Figure 3 for the Voyager 2 encounter conditions. These conditions are 
appropriate for the observations of Brown (1981) . The motion of lo in the 
torus caused by the oscillation of the plasma torus about the satellite orbit 
plane is illustrated and causes a variation of a factor about 6 in the oxygen 
lifetime at the satellite location. In the present model the centrifhgal 
plane and the satellite orbit plane are assumed coincident so that no 
oscillation of . the plasma is included. 
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To model the two-dimensional sky-plane intensity map of the 6300 X, 

0 0 

1304 A and 880 A emissions, many atoms are emitted from lo and the lifetime 

and excitation rates along their orbits are appropriately calculated and 

spatially tabulated. For model results in this paper, 1298 individual orbits 

are emitted radially and isotropically from lo's ezobase (assumed 2600 km in 

radius), each with an initial speed of 2.6 km s-^. This value of 2.6 km s-*^ 

represents a likely value for the mean speed of an initial velocity dispersion 

and was chosen because it has been the value best suited for modeling the 

spatial morphology of the bright lo sodium cloud (Smyth and HcElroy, 1978; 

Smythf 1979 and 1982) . The implicaton of velocity dispersion will be 

discussed later. Each trajectory in the model is associated with an ensemble 

of oxygen atoms emitted from a surface element of the satellite exobase. In 

model calculations, the flight time of atom trajectories must be selected 

sufficiently long to insure that a steady state description of the neutral gas 

cloud has been achieved. In all model results to follow, flight times of SOO 

hours and 1000 hours were suitable for the Voyager 1 and Voyager 2 encounter 

conditions respectively. The flux of the oxygen atoms emitted from lo is 

0 

determined by requiring that the calculated and measured 6300 A intensities 


are identical 
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3 !i Model Results 

The spatial morphologies of the emission intensities* neutral density and 

ion creaton rate for the lo atomic oxygen cloud have been calculated using the 

model described in Section 2, Because the only detection of the cloud has 

0 

been the observation of its 6300 A intensity (Brown* 1981)* model calculations 
to be presented will emphasize results at this wavelength and only briefly 

discuss intensity results for the wavelengths of 1304 X and 880 A. 

[■' 

Model results for the 6300 X intensity on the shy plane* appropriate to 

the detection observation of Brown (1981)* are shown in Figure 4 where Voyager 

2 plasma encounter conditions have been assumed. The two locations of Brown's 

observing slit are indicated* centered at S.O and S.9 B.. Bis intensity 

J 

measurement of 8 Rayleighs is an average of sis 30 minute exposures* two taken 
at the S.O slit location. This average intensity is achieved by our model 
calculation for an isotropic satellite escape flux of l.S x 10^ atoms cm*** s-^ 
or a total source rate of 6.2 x 10^* atoms s-^. The intensity contours in 
Figure 4, labeled in Rayleighs* ate evaluated for this value of the oxygen 
source. Using the hotter plasma conditions at the time of the encounter of 
Voyager 1 with Jupiter* the neutral source rate is a factor of two larger* 
i.e.* 1.2 X 10*’ atoms s-*. 

In Figure 4* the oxygon cloud is scon to completely oncirolc Jupiter, 

having a non-uniform intensity peaked near lo and also near the elongation 

point of the satellite orbit. If Brown bad used only one slit location and 

had centered it more optimally at 5.72 R. rather than at 5.9 R.* the model 

<3 J 

calculation predicts that a signal of 14.8 Rayleighs would have been measured. 

0 0 

The spatial variations of the 1304 A and 880 A intensities on the sky plane 
are similar in pattern to those of the 6300 X intensity in Figure 4. The 



14 


I 
I 

,i 

I 

0 

absolute intensity of the 1304 A omission is comparable to the 6300 
intensity for the Voyager 1 encounter plasma conditions and about two-thirds 

of this value for the Voyager 2 encounter conditions. The absolute intensity 

0 0 
of the 880 A emission is about a factor of five smaller than the 6300 A 

intensity for the Voyager 1 encounter conditions and about a factor of ten 

smaller for the Voyager 2 encounter conditions. In Table 2, the 6 Rayleigh 

0 

upper limit of the 1304 A intensity determined by long-integration-timo 
measurements of the lUE satellite UV instrument exceeds, because of its 

I ^ 

relatively large viewing aperature, the corresponding model simulated value. 

0 

The other upper limit for the 1304 A emission intensity of 25 Rayleighs in 
Table 2 was established by the Voyager 1 UVS instrument when its much smaller 
rectangular viewing slit (0.5 ^ x 0.05 
Near lo the model calculated emission brightnesses are approximately symmetric 

north, south* east and west of the satolHte inside a radius of about 0.5 Rj, 

o . , 

Model calculated values of the 13Q4 A intensity seen by the UV slit are 

14 Rayleighs for Voyager 1 encounter conditions and 6 Rayleighs for the 

Voyager 2 encounter conditions, both below the measured upper limit. The 

10 Rayleigh upper limit of 88Q X emission intensity in Table 2 is also well 

above the model calculated value. 

0 

The morphology of the 6300 A intensity viewed from above the satellite 
orbit plane and corresponding to the model calculation of Figure 4 is shown in 
Figure 5. The intensity is brightest near lo and is moderately bright ahead 

t 

of the satellite and just inside its orbit, where the temperature is still 
sufficiently hot to excite the line emission and where neutral oxygen is 
relatively abundant. The distribution of neutral oxygen is shown from the 
same viewing perspective in Figure 6. Oxygen atoms are seen to be 


Rj) was centered on the satellite. 
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concentrated primarily inside of the satellite orbit where the electrons are 
too cold for ionization of the cloud. Charge exchange not included in these 
model calculations, will, however, modify to some extent the neutral 
distribution inside of lo's circular orbit. The calculated instantaneous ion 
creation rate produced by electron Impact ionization of the oxygen cloud and 

corresponding to the results of Figure 6 is presented in Figure 7. Ahead of 

! 

lo, the ions are mostly created just inside of the satellite orbit, whereas 

I 

behind lo, most ions are created outside of the satellite orbit and somev/hat 
nearer the satellite where ionization lifetime has its minimum value (see 
Figure 3) . This overall spatial pattern is a direct consequence of the motion 
of neutrals in the gravitational force fields of lo and Jupiter and the 
spatial structure of the lifetime of oxygen in the plasma torus. 

A radial profile for the abundance of neutral oxygen and for the ion 
creation rate can be produced by integrating the individual three-dimensional 

i 

distributions over the angular dimension around Jupiter and the vertical 
dimension normal to the satellite plane. These radially averaged profiles, so 
constructed for radial integration- intervals of 0.1 Jupiter radii, are 
presented in Figure 8. The profiles were computed for the same Voyager 2 
plasma conditions and satellite emission conditions assumed in Figures 4-7. 
llie oxygen cloud contains 2.28 x 10^* atoms and the total ion source rate of 
6.20 X 10^^ ions s-^. The abundance of neutrals is seen to decrease rapidly 
outside of lo's orbit while the ion creation rate is somewhat inflated because 

t 

of the minimum lifetime of oxygen at 7.3 R. (see Figure 3). Inside of lo's 

J 

orbit the neutrals are more abundant and the ion creation rate drops rapidly 
because of the precipitous decrease in the electron temperature (see Figures 1 


and 3) 
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4. Implications for the Plasma Toms 

Tbe calculations presented in Section 3 for the atomic oxygen cloud are a 

( 

necessary first step in better understanding several physical processes that 
take place in tbe lo plasma torus. The rate at which oxygen atoms are lost by 

the neutral cloud may be directly used to provide a lower limit for tbe ion 

• » • 
loading rate N» plasma mass loading rate M and ion energy input rate E to the 

plasma torus. Tbe oxygen results determine only lower limits for these three 

quantities since other yet unidentified direct plasma sources may also exist 

for the magnetosphere as well as other neutral gas clouds. The lower limit 

for the ion loading rate N may then be used to estimate an upper limit for the 

ion residence time (i.e., an ion-diffusion loss time) in the plasma torus apd 

thereby constrain the form of the plasma diffusion coefficient for the 

magnetosphere. 

Tbe procedure to be utilized here for atomic oxygen may be generalized to 

k 

several neutral gas clouds (0, S, etc.) where charge exchange is included and 

requires that only the neutral loss rates L. and tbe neutral production rates 

J . . 

P. be calculated for each of the neutral species j. In tbe general case 
J 

N=-2(P.-L.) (3) 

j J J 

M = -Zm (P -L ) (4) 

j J J J 

E » -a/2 mjjy|. V. - Lj] (5) 

i 

wbere is the mass of the j neutral species* is the average thermal 
^ velocity of an ion previous to being converted to the j neutral species 
(which then escapes the Jupiter system)* and is the relative velocity of an 
ion with respect to tbe rotating magnetosphere upon being produced by 
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til • • 

ionization of the j neutral species. Ilie expression for N and M follow 
directly from conservation ref^uirements between the ion> electron, and neutral 
components, and the expression for E assumes that ion-ion charge exchange 
reactions do not significantly provide a net production or loss of plasma 

torus energy. In contrast to the plasma mass loading rate (4), a mass loss 

« 

rate for the neutral cloud may also be defined: 

i =Zm.L, (6) 

i c j j j 

■f 

f- 

The results of Section 3 represent the simplest case for the quantities (3) - 
^) , where the neutral, loss rate depends only upon electron impact ionization 
of atomic oxygen and where the neutral production rate is zero: 


N = L 

0 

(7) 

t • 

M *= M = m L 

(8) 

coo 

E = 1/2 m„V»L 
0 R o 

(9) 


The radial profile for L calculated for the plasma encounter conditions 

I 'O' 

of Voyager 2 is given in Figure 8 by the oxygen ion creation rate and produces 
Spatially integrated source rates for N, H, and E of 6.2 x 10^* ions s-^, 

16.6 kg s-^, and 2.7 x lO^'" watts respectively. In calculating this last 
quantity, a value of 56.8 km s-^ for was assumed and determined by the 
relative motion of lo and the rotating magnetosphere. An upper limit for the 

ion-diffusive loss time v in the plasma torus may be estimated from N as 

1 

follows 



T 


( 10 ) 
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where N is the number of ions in the spatial element in which N is calculated. 
Using the spatially integrated source rate and an estimated total of 
2.7 X 10^^ ions, the averaged ion-diffusive loss time is constrained to 

SOO days. This relatively large value is reduced significantly when 

« 

additional contributions to N, roughly estimated below, are included for a 

( 

sulfur gas cloud and for charge exchange processes, ^e spatially integrated 
value of 2.7 x 10^<* watts for E, which represents less than 1% of the 
3 X 10^> watts radiated by the UV torus, will likewise be significantly 

; I ■ 

increased by these additional contributions estimated below. For the plasma 

encounter conditions of Voyager 1, the spatially integrated source rates for 

» • » 

N, M, and E are a factor of two larger than those given above, if it is 
0 

assumed the 6300 A intensity at this earlier time had the same brightness. 

The strength of the conclusions drawn above from the results of Section 3 
are limited to some extent by several uncertainties in modeling. Inaccuracies 
in excitation cross-sections as well as in data that specify. the plasma torus 
properties could produce changes to increase or decrease the results. The 
lack of knowledge about the oxygen source characteristics of lo (i.e., the 
importance of higher velocity components in the dispersion of the emission 
velocity distribution and its exospheric distribution) could raise the overall 
satellite flux required to maintain the observed 6300 X brightness, but 
competing processes make the overall effect difficult to estimate without 
further observational information about the angular extent of the atomic 
oxygen cloud. On the other hand, inclusion of an oscillating torus in the 
model would produce a longer average lifetime for atomic oxygen and would 

0 

therfore likely reduce the satellite flux required to match a given 6300 A 
observation. The two oppositely directed effects may to first order cancel 


each other 
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Significant modifications of the conclusions drawn above from the results 
of Section 3 may, however, be anticipated from two factors not included in the 
model: (1) the likely existence of a neutral sulfur cloud for lo and (2) 

the importance of several neutral-ion and ion-ion charge exchange reactions in 

I 

the plasma torus involving the species 01, Oil, OIII, OIV, SI, SII, Sill, and 
SIV. The presence of charge exchange reactions involving slow neutrals and 
fast corotating ions of the form 

! 


neutral + 

ion 


ion + 

neutral 

(11) 

neutral 

ion 


ion + 

ion 

(12) 


both enhance and couple the oxygen and sulfur loss rates as well as provide 

substantial neutral production rates. An improved model incorporating this 

coupling is then required to accurately describe the neutral clouds and their 

impact on the magnetosphere. Lacking this coupling in the present model an^ 

the necessary information about the presence of neutral sulfur in the Jovian 

environment, this impact may be roughly estimated by scaling the above 

spatially integrated model results by the model calculation of Brown, 

Shemansky and Johnson (1983) . In their calculation, the ratio of the overall 

oxygen to sulfur loss rate was assumed to be two, which is the ratio for a 

satellite loss rate based upon the source molecule SO,. Using t^^s 

calculation as a guide, the spatially integrated value of the oxygen loss rate 

calculated above for electron impact ionization only is enhanced by a 

factor of eight by charge exchange reactions, providing a scaled value of 

. 5 X 10^7 01 atoms s-^. The spatially integrated values of the sulfur loss 

rate L is then one-half this value or 2.5 x 10^^ SI atoms s-^. These two 
s 
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loss rates are equivalent to satellite fluxes of 1.2x10^0 01 atoms cm-‘s-^ and 

6 .0 X 10” SI atoms cm-” s-” respectively. Uxe spatially integrated values of 

the oxygen production rate P and sulfur production rate P may likewise be 

o s 

scaled using the results of Brown, Shemansky and Johnson (1983) by a factor of 
approximately 0.6 and 0.3 times their respective loss rates yielding values of 

3 .0 X 10”7 01 atoms s-” and 0.75 x 10”^ SI atoms s-”. These production rates 
arise from the charge exchange reactions (11) and also from electron 
recombination reactions. 

' Using these approximations for the production and loss rates for neutral 
oxygen and sulfur in the expressions (3) through (6), more realistic values 
for the ion loading rate, the plasma mass loading rate, the cloud mass loss 
rate, and the ion energy input rate for the plasma torus may bo calculated and 
expressed as 

N=0.75(L^) (13) 

M = l.l(m^L^) (14) 

=• 2.0(m^L^) (15) 

E = 2.0(l/2 m V”)L - 0. 9(1/2 m V”)L (16) 

0 K O O i O 

where the spatially integrated value of is 5 x 10” ^ ions s-”. The effects 

of charge exchange have thus been to modify the simple dependence exhibited 

« • « • 

between N, M, M and E in the earlier derived expressions (7), (8) and (9) 
c 

which were based only on electron impact inonization. The spatially < 
integrated value of the ion loading rate (13) is then 3.85 x 10” ? ions s-” and 
' implies a value for the ion diffusive loss time (10) of 200 days, assuming a 
density of 1380 cm-” (appropriate to the Voyager 2 encounter) and a plasma 
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il 
1 1 

■j 

f 
I 

torus of radius centered on lo's orbit. The spatially integrated value of 
the plasma mass loading rate (14) is 146 kg s-^ of which S3 kg s-^ M supplied 
by oxygen atoms and S3 kg s-^ by sulfur atoms. In contrast, tl^e spatially 
integrated value for the mass loss rate of the neutral cloud (IS) is 
266 kg s-^. The spatially integrated value of the ion energy input (16) is 
3 .9 s 10^^ watts, using a value for V,p = 26.9 km/s appropriate to the average 
temperature of 60 eV reported by Brown (1982) for S^. or is 4.3 x 10*^*’ watts 
using a smaller thermal speed for the ions of 7 .6 km s-^. These two values of 
the ion energy input rate represent respectively 13% and 14% of the 
3 X 10^^ watts radiated by the plasma torus in BV wavelengths. This ion 
energy input is similar in magnitude to the lo correlated energy source 
discovered by Sandel and Broadfoot (1982) which supplies about 20% of the UV 
radiated power. 

To determine the spatial variations of the ion number, mass and energy 
injection rates in the plasma torus supplied by the neutral clouds, more , 
refined model calculations (which are currently under development) will be 
required. The overall spatial character of these processes for oxygen or 
sulfur may however be understood as being controlled through three channels if 
charge exchange reactions are considered. Electron ionization in the present 
model would tend to dominate in radial regions more distant than 6.7 Bj. In 
the radial interval S.6 to 6.7 Itj. roughly equal contributions from electron 
impact ionization and charge exchange are expected. Charge exchange reactions 
will dominate inside of 5.6 Bj where the electron temperature decreases 
rapi41y. Charge exchange processes are therefore expected to modify the 
' radial distribution of 01 atoms and the ion creation rate shown in Figure 8 
which are calculated only for electron impact ionization. For radial 
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displacements larger than 5.9 Uj, the rather sharp reduction in the 01 
population and the substantial ion creation rate out to 7.5 Bj. are caused by 
the rapidly rising elecron ionization rate coefficient and these features will 
have the same character when charge exchange is included. The drop in the ion 
creation rate through two decades in the inward direction from 5 .9 to 
5.5 B^ occurs because of the precipitous drop in the electron temperature in 
that region. The inclusion of charge exchange reactions will provide a now 
ion production mechanism inward of 5.9 Bj. which will enhance the ion creation 
rate and will also sharpen the neutral atom distribution by reducing the 
relative number of atoms in the 4.5 to 5.5 B^ region. The extent of the 
oxygen cloud will also be confined more closely to the position of lo and the 
radial distribution of oxygen may then tend to show a slightly higher density 
peak than that shown in Figure 8, because the Brown (1981) observation was 
made ~90^ downstream from the position of lo. 

The above estimates for the ion loading rate, the mass loading rate and 
the ion energy loading rate for the plasma torus have as their source the 
neutral gas clouds of Ip. Although it is almost certain that lo is the 
ultimate source of these ions, other sources of neutrals or direct sources of 
ions that have yet to be discovered might also exist in addition to the ions 
created by the neutral gas clouds. If the other sources of ions are 
significant, the above estimates of N. M and E are only lower limits. It is 
difficult to assess the possible magnitude of direct ion sources since 
knowledge of lo's atmosphere and the mechanisms which deliver the material to 
the torus is very limited. Only oxygen ionization processes that occur in the 
' near vicinity of lo have been limited by direct observation and predictable 
emission efficiencies (Shemansky 1980). The Shemansky (1980) upper limit on 
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direct ion injection from lo is roughly 10^ ^ based on ion detection of the 

0 0 
01 1304 A line in Voyger 1 observations. Attempts to observe the 1304 A line 

by Bfoos and Clarke (1981) using the lUE satellite produce about the same 

value* 10*7 8-*, as an upper limit for the Voyager 1 encounter epoch. It is 

likely that the source molecule supplying the neutral gas clouds is SO, (Pearl 

et. 1979j Johnson et al . . 1979) so that direct escape of SO,* SO,^ or its 

chemical fragments (SO* 0,* etc.) may occur and not have as yet been detected. 

The assessment of the importance of these molecules and ions as a source of 

plasma for the magnetosphere will* however* require that additional 

inforaation be acquired. 
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5 . Concluding Remarks 

A model bus been presented for the neutral oxygen cloud of lo based upon 

atom-electron impact excitation and ionization processes in the plasma torus. 

The flux of oxygen atoms emitted fi'om lo was determined by requiring that the 

measured (Brown 1981) and the calculated 6300 A intensities be identical. By 

this meams satellite emission fluxes of 3.0 x 10* cm-* s-* and 

l.S X 10* cm-* S-* were required for plasma conditions appropriate to the time 

0 

of encounter of Voyager 1 and Voyager 2 with Jupiter. The 6300 A observations 
of Brown (1981) wore more appropriate to the Voyager 2 plasma conditions for 
which spatially integrated ion loading rates, plasma mass loading rates and 
ion energy input rates were calculated to be 6.2 x 10** ions s-*, 16.6 kg s-*. 
and 2,7 x 10** watts respectively. This value of the ion loading rate implies 
an ion-diffusive loss time of SOO days. 

Enhanced values of the satellite emission flux, ion loading rate, mass 
loading rate and ion energy input rate resulting from the presence of a likely 
existing sulfur cloud and from the impact of charge exchange reactions were 
then roughly estimated based on the results of the oxygen model and the charge 
exchange calculations of Brown. Shemansky and Johnson (1983). The upward 
adjusted value of the satellite flux for atomic oxygen was 1.2 x 10** cm-* s-* 
and an additional satellite flux of 6.0 x 10* cm-* s-* was assumed for the 
likely existing sulfur cloud of lo. This represents a neutral mass loss for 

the satellite of ~270 kg s-*. The corresponding adjusted values for the 

/ 

spatially integrated ion loading rate, plasma mass loading rate, and ion 
energy input rate were then 3 .8 x 10** ions s-*. ~1S0 kg s-* and 
(3.9 - 4.3) X 10** watts respectively. This value for the ion loading rate 
implies an ion-diffusive loss time of ~200 days. The value of the ion energy 



input rate represents about 13 - 14% of the 3 x 10^* watts radiated in the UV 
by the plasma torus and is comparable to the 20% energy input level associated 
with the lo correlated energy source discovered by Sandel and Broadfoot 
(1982) . The remaining 80% of the energy input has been associated by 
Shemansky and Sandel (1982) with a local time asymmetry in the electron 
temperature. 

Inclusion of charge exchange reactions in future models of the atomic 
oxygen cloud is important both in refining the values of physical quantities 
roughly estimated above and in providing an interpretative base for new 
observations of the neutral and ion species of the magnetosphere. Such a 
model is presently under development. New observations of the angular 
distribution of atomic oxygen around Jupiter would be particularly useful in 
refining the value of the oxygen flux emitted by lo. Observations to detect 
an atomic sulfur cloud are likewise desirable in order to establish a 
satellite emission flux for sulfur and thereby better understand processes in 
the local atmosphere of Io» 
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Table 2 


Electron Impact Ionisation R^te and Electron Impact 
Excitation Rates at 6300 k, 1304 A and 880 A for Atomic Oxygen 


Electron 

Temperature 

(eV) 


Ionization 
Rate^ 
(cm* sec-*') 


6300 A 
Excitation 
Rate® 
(cm® sec-*) 


1304 A 
on 

(cm* sec-*) 


Excitati 

Rate® 


880 A 


Excitation 
Rate® 
(cm* sec-*) 


o.s 

6.80 

X 

10-** 

3.44 

X 

10-** 

1.99 

X 

io-*» 

6.23 

X 

10-** 

T.O 

7.86 

X 

10-** 

3.20 

X 

10-*» 

3.61 

X 

10-** 

1.15 

X 

10-** 

2.0 

1.02 

X 

10-** 

1.09 

X 

10- » 

5.50 

X 

10-** 

1.79 

X 

10-*» 

S.O 

1.10 

X 

10- » 

2.40 

X 

10- » 

1.27 

X 

10- * 

1.71 

X 

10-*o 

10.0 

7.18 

X 

10- « 

2.92 

X 

10- » 

3.01 

X 

10- » 

8.45 

X 

10-*° 

20.0 

2.21 

X 

10-8 

2.72 

X 

10- » 

6.68 

X 

10- * 

1.94 

X 

,10-» 

50.0 

4.93 

X 

10-» 

1.70 

X 

10- * 

9.04 

X 

10- » 

3.15 

X 

10- * 

100.0 

6.63 

X 

10-* 

9.25 

X 

10-*o 

9.45 

X 

10- * 

3.55 

X 

io-» 

200.0 

7.54 

X 

10-* 

4.25 

X 

10-*» 

9.00 

X 

10- » 

3.55 

X 

10- » 

500.0 

7.59 

X 

10-* 

1.28 

X 

10-*« 

7.70 

X 

10- » 

3.16 

X 

io-» 

1000.0 

7.03 

X 

10-* 

4.84 

X 

10-** 

6.54 

X 

10- * 

2.74 

X 

10- » 

2000.0 

6.19 

X 

10-* 

1.77 

X 

10-** 

5.40 

X 

10- » 

2.29 

X 

10- » 

5000.0 

4.97 

X 

10-* 

4.56 

X 

10-*» 

4.07 

X 

10- * 

1.75 

X 

10- » 


a. Based on ionization cross-section of Brook., Harrison and Smith (1978). 

b. Based on 01 -t- e calculations by Shemansky (unpublished); collision 

0 

strength data for 6300 A emission were obtained from calculations of 

< 

Henry, Burke and Sinfailam (1969), measurements by Stone and Zipf (1974) 

0 

(see Shemansky (1980) for scaling factor) for 1304 A, and theoretical 

0 

estimates for 880 A. 
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Figure Captions 

Figure 1. Electron Temperature Radial Profiles for the lo Plasma Torus,' The 
two electron temperature profiles adopted for model calculations 
are discussed in the text. The plotted points are i.n situ 
measurements from Scudder, Sittler and Bridge (1981) , 

Figure 2. Electron Number Density Profiles for the lo Plasma Torus. The 

different number density profiles in the centrifugal plane of the 

I 

plasma torus are discussed in the text. The two profiles adopted 
for model calculations are indicated by dashed lines for radial 
displacements less than 5.8 R. and by solid linos for larger radial 
displacements. 

Figure 3. Electron Impact Ionization Lifetime (hours) for Atomic Oxygen in 

the lo Plasma Torus. The lifetime is calculated using the plasma 

conditions adopted for the Voyager 2 encounter with Jupiter. 

Ionization cross-sections are based upon the experimental ‘ 

measurements of Brook. Harrison and Smith (1978) . 

0 

Figure 4. lo Oxygen Torus 6300 A pnissipn Intensity, The satellite location 

and the calculated intensity contours are appropriate to the 

detection observation of Brown (1981) . Voyager 2 encounter plasma 

conditions and a satellite flux of 1,5 x 10* oxygen atoms cm-* s-* 

were assumed. See text for discussion. 

0 

Figure S, 6300 A Oxygen Emission Intensity. The intensity morphology of the 

I 

lo oxygen cloud viewed from above the satellite orbit plane is 
calculated for the same plasma and satellite emission conditions of 
Figure 4. The intensity contours have values, in units of 
Rayleighs from the outer to inner contour, of 0.045. 0.32. 0.59 and 
1.13 respectively. 
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Figure 6 


Figure 7 


Figure 8 


lo Atomic Oxygen Torus. The column density morphology viewed from 

above the satellite orbit plane is calculated for the same plasma 

and satellite emission conditions of Figure 4. The column density 

contours have values, in units of cm* from the outer to the inner 

contour of 6 ,6 s 10*, 6 .7 x 10*-®, 1.3 x 10*'*^ , 1.9 x lO*-*- and 

3 .1 X 10^^ respectively. 

A. 

. ' lo Oxygen Ion Creat^n Rate. The ion creation rate morphology 
viewed from above the satellite orbit plane is calculated for the 
same plasma and satellite emission conditions of Figure 4. The ion 
creation rate contours have values, in units of ions pm-* s-* from 
the outer to the inner contour, of 2.8 x 10^, 2.0 x 10*, 5.3 x 10*, 
8.6 X 10* and 1.4 x 10* respectively. 

Radial Profiles for Oxygen Abundance and Ion Creation Rate, The 
profiles are calculated for the same plasma and satellite emission 
conditions of Figure 4. See text for discussion. 






Electron Number Density (cm*^) 
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Figure 2 



Vertical Displacement from the Centrifugal Equator (Rj) 
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01 Electron Impact Ionization Lifetime in the lo Plasma Torus 

(Voyager 2 Plasma Conditions) 



Radial Displacement from Jupiter (Rj) 


Figure 3 
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